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, Abstract. Our Galaxy harbours a large population of X-ray sources of intermediate to low 

• . X-ray luminosity (typically Lx from 10"^ to 10'"* ergs"'). At energies below 2keV, active 

i-J^ ■ coronae completely dominate the X-ray landscape. However, the nature and the properties 

of Galactic sources detected at energies ^ 2 keV is much less constrained. Optical follow-up 
JL , spectroscopic observations show that in addition to cataclysmic variables (CVs) and very 

; I ■ active stellar coronae, massive stars (colliding wind binaries, quiescent high-mass X-ray bi- 

■ naries and y-Cas analogs) account for a sizable fraction of the Galactic hard X-ray sources 
I at medium flux (Fx ^ 10"'^^ ergcm"^ s"'). Cross-correlations of the 2XMM-DR3 catalogue 

I I. with 2MASS and GLIMPSE confirm the presence above 2keV of a large population of 

coronally active binaries, probably of the BY Dra and RS CVn types, in addition to many 
' distant and absorbed massive stars. We also report the results of a specific optical identifi- 

^ , cation campaign aimed at studying the nature of the optically faint hard X-ray sources and 

04 ■ at constraining the surface density of black holes (BHs), either isolated and accreting from 

the interstellar medium or in quiescent binaries. Not astonishingly, most of our sample of 
14 optically faint and X-ray hard sources are identified with CVs and Me stars. We do not 
, find any likely counterpart in only three cases. Our observations also allow us to put an 

■ upper limit of 0.2 BH deg"^ at Fx= 1.3 10"'^ ergcm"^ s"' in directions toward the centre of 

■ our Galaxy. This implies a combined Bondi-Hoyle and M to Lx efficiency of accretion onto 
I black holes of less than 10"'. 

• • ■ Key words. X-ray: stars - X-ray: black holes - X-ray: populations 
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5h I 1 . Introduction ergy emission and cover a wide range of spec- 

^ ■ tral properties and X-ray luminosities. Active 

Our Galaxy presents a rich and contrasted ^j^jj^^ ^^^^^^^ ^^^^ numerous 

X-ray landscape. Many different classes of soft (kT ~ 0.5 keV) X-ray sources encountered 

Galactic objects can be the source of high en- q^^^^^-^ latitudes. Their X-ray lumi- 



nosities are in the range of 10^^ to 10^' erg s ' 
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With Lx typically higher than 10^^ ergs"', 
classical High and Low-Mass X-ray binaries 
(HMXBs, LMXBs) occupy the bright end 
of the Galactic X-ray luminosity function. 
However, X-ray surveys have also shown the 
presence of a large population of relatively 
hard Galactic sources with Lx in the range 
of 10-" to a few 10-''* ergs"'. The nature and 
overall properties of these intermediate X-ray 
luminosity sources remains badly understood. 
Because of their faintness, these populations 
can only be studied in our Galaxy. Although 
CVs and coronally active binaries are likely 
to account for a large part of this population, 
many other kinds of objects have also been 
identified in this range of Lx, namely X-ray 
binary transients in quiescent states, magnetic 
OB stars or colliding wind massive binaries. 
Among the most interesting objects that could 
appear in this Lx regime are isolated compact 
remnants of early stellar formation accreting 
from the interstellar medium and low luminos- 
ity stages predicted by binary evolution models 
(e.g.. Be + white dwarf, wind accreting neutron 
star binaries or precursors of LMXBs). 

Low-luminosity hard X-ray sources are 
also hkely to be important contributors to the 
Galactic ridge X-ray emission (GXRE; see 
e.g. Worrall et al. 1982). The nature of this 
extended and apparently diffuse emission re- 
mains debated. Its X-ray spectrum displays a 
prominent emission hne at 6.7 keV and resem- 
bles that of a thin thermal plasma with tem- 
peratures of 5-lOkeV (Koyama et al. 1986). 
However, such a hot interstellar medium can- 
not remain bound in the Galactic gravitational 
potential well and its presence would require 
that unrealistic ally powerful sources of hot 
plasma concur to replenish it continuously. 
The close similarity of the spatial distribu- 
tion of the GXRE with near-infrared emission 
(Revnivtsev et al. 2006) also favours an ex- 
planation in terms of unresolved emission of 
many low-Lx sources. A deep Chandra ob- 
servation of one of the brightest ridge emis- 
sion close to Galactic Centre resolved at least 
80% of the diffuse emission into point sources 
at energies 6-7 keV (Revnivtsev et al. 2009). 
However, Ebisawa et al. (2005) reach opposite 
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conclusions based on another deep Chandra 
observation at / ~ 28.5°. 

The goal of this paper is to explore the na- 
ture of low to intermediate X-ray luminosity 
sources encountered in the Galaxy and shining 
in the hard (2-12keV) range. Our work is based 
on spectroscopic identifications obtained at the 
telescope and identifications derived from the 
cross-correlation of XMM-Newton serendipi- 
tous sources with large optical and infra-red 
archival catalogues. 



2. The X-ray content of the 
XMM-Newton Galactic Plane 
Survey 

The Survey Science Centre (SSC) of the 
XMM-Newton satellite has recently reported 
results from an optical campaign aiming at the 
identification of the brightest X-ray sources in 
the XGPS (Motch et al. 2010). The ~3deg2 
area surveyed is located at Z = 20°, ^ = 0° 
(Hands et al. 2004). Among the 29 hard (2- 
lOkeV; S/N > 3) sources investigated optically, 
six are identified with massive stars possibly 
containing an accreting component or being 
powered by colliding winds, three are identi- 
fied with CVs, two with low-mass X-ray bi- 
nary candidates and six with stars. At Fx ^ 
10"'^ ergcm"^ s"' (2-lOkeV), a large fraction 
of the expected Galactic source population is 
positively identified. Active coronae account 
for ~ 10% of the expected number of Galactic 
sources in the hard band. 



3. 2MASS and GLIMPSE 
Identifications 

The advent of high quality photometric infra- 
red surveys covering large areas offers a unique 
opportunity to statistically identify and charac- 
terise hard Galactic X-ray sources. Using the 
method outlined in Pineau et al. (2011), we 
computed probabilities of identification with 
the 2MASS and GLIMPSE catalogues for all 
2XMM-DR3 sources located within 3° from 
the Galactic plane. About one quarter of the 
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-38,000 2XMM-DR3' low b entries have a 
2MASS match with an individual probability 
higher than 90%. At this level, the expected 
number of spurious matches is ~ 1.3% (Motch 
et al. 2010). Figs land 2 show the distribu- 
tion in EPIC pn hardness ratios of the 2XMM 
DR3 sources with and without 2MASS coun- 
terparts. Sources matching 2MASS entries are 
clearly much softer in X-rays. The peak in the 
HR2 histogram is consistent with optically thin 
thermal emission with kT ~ 0.5 keV undergo- 
ing an absorption with logNe ~21.5 while the 
HR3 distribution indicates the presence of a 
harder X-ray component (kT~lkeV). Most 
2MASS identifications are thus likely active 
coronae. The few soft XMM-Newton sources 
without 2MASS entries are likely Me stars be- 
ing too faint in the infrared to be listed in the 
2MASS catalogue. 

We show in Figs. 3 and 4 the distribution 
of EPIC pn HRs with the H-K colour. Since 
the intrinsic stellar H-K colour index remains 
within -0. 1 to +0. 1 from O to K spectral types 
and all luminosity classes (see e.g. Covey et al. 
2007), the H-K colour mainly reflects interstel- 
lar absorption. No single X-ray energy distri- 
bution can account for the overall HR/E(H-K) 
relation. The most absorbed sources appear to 
be also the intrinsically hardest ones. The hot- 
ter X-ray temperature of young stellar coronae 
and active binaries combined with their higher 
luminosity make them detectable up to larger 
distances than older and X-ray softer stars. In 
particular, the bulk of the stars well detected 
above 2 keV (i.e. appearing in Fig. 4) have HR3 
consistent with active binaries of the BY Dra or 
RS CVn type. 

However, many very hard X-ray sources 
have high probability 2MASS identifications. 
Their number is significantly larger than ex- 
pected from spurious matches. Assuming that 
their red H-K is due to interstellar absorption 
(Nh ~ 10^2 for E(H-K) = 0.3) yield dis- 

tances larger than 3kpc (for a mean particle 
density of n=l). Sources located above the 
F = 2 powerlaw curve in the H-K versus HR3 




-0.5 0.5 

pn HRS (0.5-1.0 keV)/{ 1.0-2.0 keV) 



Fig. 1. Distribution of EPIC pn hardness ratio HR2 
= ([1.0-2.0] -[0.5-1.0])/[0.5~2.0keV]) for \b\ < 3° 
2XMM-DR3 sources with err(HR2) < 0.1. Black: 
sources having a > 90% probability to be associ- 
ated with a 2MASS entry; red: sources without any 
2MASS entry within a combined X-ray -I- 2MASS 
3cr error radius. 
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' The third release of the 2XMM catalogue 
was published in April 2010. see http://xmmssc- 
w w w. s tar. le . ac . uk/C atalogue/2XMMi-DR3/ 



Fig. 2. Same as Fig.l for EPIC pn HR3 = ([2.0- 
4.5]-[1.0-2.0])/[1.0-4.5keV]. 



diagram shown in Fig. 4 have a H magnitude of 
about 1 1 and a H-K around 0.45. At a distance 
of 3kpc, the absolute H magnitude is <-1.7 
suggesting that these 2MASS identifications 
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Fig. 3. H-K colour versus EPIC pn HR2 for all 
\b\ < 3° 2XMM-DR3 sources with err(HR2) < 0.25, 
err(H-K)<0.1 and probability of 2MASS identi- 
fication >90%. The first three lower dark (red) 
lines show the expected E(H-K) versus HR2 rela- 
tion assuming 2-T thermal emission from a 1.9 Gyr, 
300 Myr and 30Myr active stars (Guedel et al. 
1997). The two hardest dark (red) relations corre- 
spond to AY Get, a typical BY Dra binary (Dempsey 
et al. 1997) and to the RS GVn star WW Dra 
(Dempsey et al. 1993). The two light (green) lines 
correspond to power laws with photon indices T of 
and 2. The RS GVn relation is almost superposed 
on the r = 2 line. Positions of the WR star XGPS- 
14 (Motch et al. 2010), (blue stars) and of a number 
of INTEGRAL HMXBs (big magenta squares) are 
also shown for comparison. Multiple detections are 
plotted. 



could well be massive stars. Their positions in 
the HR/H-K diagram are indeed comparable to 
those of HMXBs discovered by INTEGRAL 
and are akin to the Wolf-Rayet (WR) XGPS- 
14 (Motch et al. 2010). Therefore, based on 
the optical and infrared spectroscopic identi- 
fication work reported by Motch et al. (2010) 
and by Anderson et al. (2011), we expect that 
many of these hard X-ray 2MASS identifica- 
tions are low-Lx HMXBs, massive stars (ei- 
ther single or in wind colliding binaries) and 
y-Cas analogs. Because of the high density 
of Spitzer sources in the Galactic plane, the 
probability to find at random a relatively bright 



Fig. 4. Same as Fig. 3 for EPIG pn HR3. WRs -l- 
INTEGRAL HMXBs are even better separated from 
coronal sources in this diagram. 

GLIMPSE source in the XMM error circle is 
rather large. Only 222 2XMM-DR3 sources, 
for which boresight correction could be carried 
out (see Watson et al. 2009), have a match- 
ing probability higher than 90%. To this we 
add 103 sources with P> 90% resulting from 
the crosscorrelation with XMM observations 
for which no boresight correction was possi- 
ble. These sources have on average larger error 
circles (Watson et al. 2009). XMM-GLIMPSE 
sources with a high probability of identifica- 
tion have a distribution in hardness ratios simi- 
lar to that seen for 2MASS. Many of the bright 
XMM-GLIMPSE sources do have identifica- 
tions with known optically bright active coro- 
nae. Among the hardest sources we find sev- 
eral catalogued WR stars, HMXBs discovered 
by INTEGRAL and a number of unidentified 
sources sharing very similar X-ray and infrared 
properties. 

4. Hard X-ray sources with faint 
optical counterparts: - A Search for 
low Lx accreting black holes 

Optically faint hard X-ray sources are much 
less constrained than those associated with lu- 
minous objects such as the massive stars dis- 
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cussed in the two previous sections. Deep in- 
frared observations have shown indeed that the 
majority of the hard X-ray sources detected 
in the Galactic Centre region are not associ- 
ated with massive stars (Laycock et al. 2005). 
Optical follow-up of ChaMPlane sources in the 
Galactic bulge (Koenig et al. 2008) suggests 
the presence of a large population of active bi- 
naries and young stellar objects with a small 
CV contribution. Our ESO-VLT optical obser- 
vations of the brightest hard X-ray sources in a 
region located at / = 20° = 0° (Motch et 
al. 2010) led to the identification of five CVs or 
LMXB candidates and of a few active coronae. 
This distribution is globally consistent with the 
locally determined X-ray luminosity function 
of faint point sources reported by Sazonov et 
al. (2006). However, we expect the contribu- 
tion of the various kinds of hard X-ray emitters 
to strongly vary with Galactic position. 

In addition, apart from toward the very cen- 
tral regions of the Galaxy, a large fraction of 
the hard X-ray sources detected at low Galactic 
latitude are background AGN. This extragalac- 
tic "contamination" depends sensitively on the 
direction of observation (Motch 2006; Hong et 
al. 2009). A further difficulty arises from the 
fact that CV X-ray spectra resemble those of 
mildly absorbed AGN and in general cannot be 
efficiently preselected on the basis of hardness 
ratios only. 

In order to increase our chances to se- 
lect genuine Galactic X-ray sources for op- 
tical follow-up, we used the signature of the 
large photoelectric absorption imprinted on 
background AGN to achieve a high rejection 
rate for extragalactic sources. Our goals were 
twofold. First, investigate the nature of this op- 
tically faint hard X-ray Galactic population and 
second, constrain the surface density of low X- 
ray luminosity black holes (BH) in quiescent 
binaries, or being isolated and accreting from 
the interstellar medium. 

Apart from three cases of long micro- 
lensing events possibly due to isolated black 
holes (Bennett et al. 2002; Nucita et al. 2006) 
and perhaps a couple of massive unseen com- 
panions in X-ray quiet binaries, all established 
or candidate stellar mass black holes (~ 60) are 
in accreting binaries (see Ziolkowski 2010, for 



a recent census). However, the actual number 
of isolated black holes (IBHs) in the Galaxy 
could be of the order of 10^ (Samland 1998; 
Sartore & Treves 2010) or even higher. A frac- 
tion of them might accrete matter from the in- 
terstellar medium at a rate high enough to be- 
come detectable in the radio, optical or X-ray 
domains. IBH X-ray properties have been in- 
vestigated in details by Agol & Kamionkowski 
(2002). Their detectabiUty depends sensitively 
on three parameters. The first one is the as- 
sumed spatial velocity distribution, which re- 
flects the amplitude of the kicks received at 
birth, and is critical for Bondi-Hoyle accre- 
tion. Jonker & Nelemans (2004) have shown 
that BH X-ray binaries display scale heights 
comparable to those of neutron star binaries. 
This suggests that at birth BHs receive ve- 
locity kicks similar to those of neutron stars. 
The second key parameter is the efficiency of 
Bondi-Hoyle accretion, which can be quite sig- 
nificantly diminished in presence of magnetic 
fields. The accretion rate can be expressed as: 



M = .1 



(1) 



with Vrei the relative velocity of the accret- 
ing object with respect to the ISM having 
sound speed c, and mass density p. Pema et 
al. (2003) argue that the value of the dimen- 
sionless parameter A which determines the ac- 
tual efficiency of Bondi-Hoyle accretion, and 
is usually assumed to be of the order of unity, 
could be as low as 10"^ or even less. The last 
important parameter is the efficiency e with 
which X-rays are generated in the accretion 
flow (Lx = eMc^). The general agreement is 
ffiat € decreases sensitively at low M due to the 
onset of radiatively inefficient accretion flows 
and as a result of an increasing fraction of ac- 
cretion energy being transformed into kinetic 
energy of possible jets. 

BH quiescent binaries display powerlaw X- 
ray spectra with photon index F between 0.9 to 
2.3 and X-ray luminosities in the range of 10^*^ 
to lO^^ergs"! (Kong et al. 2002; Hameury 
et al. 2003). Here we assume that IBHs ex- 
hibit similar X-ray spectra. As a best com- 
promise between maximum distance range and 
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survey area, we selected Galactic directions {\b\ 
< 3°) toward which the total Galactic inter- 
stellar column is higher than lO^^cm"^. For a 
typical ISM density of n ~ 1 cm"^, this Nh is 
reached at distances of ~ 3 kpc. At the time 
we started this project we relied on obser- 
vations contained in the IXMM catalogue to 
which we added several pointings extracted 
from the XMM-Newton archive. In order to 
obtain consistent HR definitions, we only con- 
sidered the EPIC pn camera. The total area ful- 
filling these conditions in our source selection 
was 20 deg2 at Fx > 5 10"''* erg cm^^ ^-i (o.2- 
12 keV), taking into account overlapping expo- 
sures. This flux corresponds to a limiting X- 
ray luminosity of 6 10""' ergs"' at Ikpc. We 
then considered all sources displaying hard- 
ness ratios consistent within the errors with 
r = 0.9-2.3 powerlaw energy distributions ab- 
sorbed by Nh < lO^^cm"^. Among the latter, 
we selected for optical follow-up 14 of the 
X-ray brightest sources having optical candi- 
dates fainter than B = 18^ and being observ- 
able from ESO La Silla during Chilean win- 
ter time. All selected X-ray sources were vi- 
sually checked and we discarded areas of high 
diffuse X-ray and optical emission as well as 
star forming regions. The set of hardness ratios 
used in the IXMM provides less energy res- 
olution in the soft bands than those used sub- 
sequently. Nevertheless, the resulting IXMM 
based source selection was found to be consis- 
tent with that based on 2XMM HRs as shown 
in Fig.5. Source selection and validation have 
been done using the XCat-DB^^ (Motch et al. 
2009), the official SSC interface to XMM cat- 
alogues developed in Strasbourg. 

Optical observations were carried out with 
the ESO-NTT + EMMI instrument from July 
31 to August 2, 2005. All spectra were ob- 
tained with Grism # 5 through a 1.5" slit. This 
setup provided a spectral resolution of 1000 in 
the wavelength range of 3800 to 7000A. 

We list in Tab. 1 a summary of the results of 
our optical identification work. Among these 



^ A fraction of the accretion luminosity should be 
emitted in the optical and infrared domain via syn- 
chrotron emission of the hot accreting plasma 

' http://xcatdb.u-strasbg.fr/ 
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Fig. 5. Distribution of candidate accreting black 
holes in the EPIC pn HR2/HR3 diagram (as defined 
in the 2XMM DR3 catalogue). The red lines de- 
fine the area containing sources with X-ray energy 
distributions expected from accreting BHs undergo- 
ing Nh < 10^-cm"-. Candidates extracted from the 
2XMM DR3 are shown as black dots. Sources se- 
lected for optical follow-up observations are plotted 
in blue together with their HRs errors. 

14 sources we identified four cataclysmic vari- 
ables, four Me stars and three active stars. 

Only three XMM sources remain uniden- 
tified with candidate counterparts fainter than 
typically W -22. We note however, that our op- 
tical spectroscopic limit is not yet deep enough 
to rule out an identification with a high Fx/F^pt 
CV such as for instance 2XMM J 18325 1.4- 
100106 (Fx = 8.5 10"'^ erg cm^^ s"', V = 23.2; 
Motch et al. 2010). From this optical cam- 
paign, we derive an upper limit of 0.2 
BH deg-2 at Fx = 1.3 10"'^ erg cm^^ §-' (0.2- 
12keV). This is ^ 20 times the IBH surface 
density predicted by Agol & Kamionkowski 
(2002) for central regions of the Galaxy such as 
that covered by our optical targets (b = -1° ,+2° 
and l--49°,+33°). According to their model 
we expect 0.17 IBH in the total area surveyed 
at Fx = 1.3 10"'^ ergcm"^ s"' and less than 3 
at a 10 times lower flux limit. It should also be 
stressed that our limiting Nh de facto excludes 
from our sample all BH located at distances 
greater than a few kpc. This negative result 
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Table 1. Summary of optical identifications. AC - Active coronae earlier than M. 



Source name X-ray flux " V mag Identification 



2XMM J135859.3-60I518 


1.87 


> 22 


UNID 


2XMM J154305.5-522709 


7.93 


20.8 


CV 


2XMM J172803.0-350039 


10.6 


18.5 


Me 


2XMM J174504.2-283552 


0.58 


17.5 


AC 


2XMM J1745 12.9-290931 


0.40 


20.9 


Me 


2XMM J175520.5-261433 


2.07 


> 22 


UNID 


2XMM J 180235. 9-23 1332 


2.87 


20.0 


CV 


2XMM J 180243.0-224 105 


1.43 


20.1 


Me 


2XMM J180913.0-190535 


2.39 


22.4 


CV 


2XMM J181003.2-212336 


8.24 


16.2 


AC 


2XMM J181857.9-160208 


1.06 


>22 


UNID 


2XMM J182703.7-1 13713 


0.65 


17.3 


AC 


2XMMJ183228. 1-102709 


0.50 


20.3 


Me 


2XMM J185233.2-I-000638 


0.62 


20.4 


CV 



" In units of 10"" erg cm'^ s"' (0.2-12 keV). Average of all de- 
tections. 



clearly illustrates the difficulty of searching 
for Galactic IBHs in the X-ray domain where 
the background of AGN and the foreground 
of hard coronal emitters and CVs completely 
dominate the population. Using the population 
parameters of Ago! & Kamionkowski (2002) 
our constraints can be expressed as A x e ^ 
210''^ /Ng with Ng the number of Galactic 
BHs in units of 10^. For a reasonable value 
of A^9 = 0.2, we derive an upper limit of « 10"^ 
on the global (Bondi-Hoyle times X-ray) effi- 
ciency. 

5. Conclusions 

In this paper, we report on several efforts 
made to characterise the nature of the serendip- 
itous XMM-Newton sources discovered in 
the Galactic plane. Dedicated optical follow- 
up observations and cross-correlation with 
archival catalogues reveal a significant number 
of hard X-ray emitting massive stars with lumi- 
nosities in the range of ~ 10^^ to ~ 10^"* erg s"' . 
This population mainly consists of Wolf-Rayet 
stars, wind colliding binaries, HMXBs in qui- 
escence and y-Cas analogs. Young stars and 
active corona binaries of the BY Dra and 



RS CVn types also contribute significantly to 
the population of sources detected at energies 
above 2keV. The optically faintest Galactic 
hard X-ray sources are mostly identified with 
cataclysmic variables. 

The typical XMM-Newton sensitivity al- 
lows us to constrain the nature of the hard X- 
ray sources of low- to intermediate- Lx up to 
distances of a few kpc. It is therefore unclear 
whether these studies can be used to assess 
the nature of the unresolved population respon- 
sible for the Galactic Ridge X-ray emission 
mostly seen at \l\ ^ 50°. 

Finally, we report on a first dedicated 
search for low-Lx black holes in quiescent bi- 
naries or in isolation and accreting from the 
interstellar medium. We derive an upper limit 
of 0.2 BH deg-2 at Fx = 1.3 10"'^^ erg cm^^ 
(0.2-12keV) in directions of the central parts 
of the Galaxy. This observational limit implies 
that the efficiency e with which X-rays are gen- 
erated in the accretion flow is $ 10"^ if one 
assumes nominal Bondi-Hoyle mass accretion 
rates. 
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DISCUSSION 

MARAT GILFANOV: Is there a similar on- 
going elfort for the soft band ? 

CHRISTIAN MOTCH: The Survey Science 
Centre of the XMM-Newton satellite has un- 
dertaken a wide area optical follow-up pro- 
gramme to identify serendipitous EPIC sources 
in different directions of the Galaxy. Most soft 
X-ray sources are identified with active coro- 
nae. The analysis of their X-ray and optical 
properties is one of the main scope of this iden- 
tification programme. 

TAKESHI GO TSURU: The X-ray spectrum 
of the Galactic ridge emission is characterised 
by iron emission lines. Do you think if the col- 
lection of point sources can explain the ridge 
spectrum including iron emission Unes ? 

CHRISTIAN MOTCH: According to 
(Revnivtsev et al. 2009), more than 80% of the 
diffuse flux at energies of 6-7 keV is resolved 
into discrete sources at / = 0.08°, b= -1.42°. 
Therefore, it seems that the combination 
of faint stellar and CV sources responsible 
for the GRXE has intrinsic iron line flux 
strong enough to account for that seen in the 
unresolved emission. 



